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F? 1 ABSTRACT 

> : 

^ ! We have derived a primordial helium abundance of Y p = 0.2477±0.0029, based 

on new atomic physics computations of the recombination coefficients of He I and 
of the collisional excitation of the H I Balmer lines together with observations 
and photoionization models of metal-poor extragalactic H II regions. The new 
atomic data increase our previous determination of Y p by 0.0086, a very significant 
amount. By combining our Y p result with the predictions made by the standard 
Big Bang nucleosynthesis model, we find a baryon-to-photon ratio, rj, in excellent 
agreement both with the rj value derived by the primordial deuterium abundance 
value observed in damped Lyman-a systems and with the one obtained from the 
WMAP observations. 



Subject headings: early universe — galaxies: abundances — galaxies: individual 
(SBS 0335-052, I Zw 18, Haro 29) — galaxies: ISM H II regions — ISM: 
abundances 
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Introduction 



This is the third paper of a se ries on the determina tion of the primordial helium abun- 
dance by unit mass, Y p . In Paper I ( IPeimbert et al.ll2002l) we studied the effect of temperature 
variations on the determination of Y p , and in Paper II (ILuridiana et al.ll2003l ) we studied the 



effect of collisional excitation of the Balmer lines in the determination of Y p . 

The determination of Y p is important for the study of cosmology and the chemical evolu- 
tion of galaxies. In particular, Y p can be used to test the standard Big Bang nucleosynthesis 
(SBBN) scenario by setting strong constraints on: a) the number of neutrino families, N u ; 
b) the variation of the neutron lifetime and of the neutron-proton mass difference, these 
constraints can be translated into constraints on the time variation of the Yukawa couplings 
and the fine structure constant; c) the variation of the constant of gravity, G; d) the presence 
of vacuum energy during BBN; and e) the presence of decayin g particles during BBN, which 



could have affected the production of the light elements (e.g. iCyburt et ajj [2005| ; ICoc et al. 
20061 . and references ther ein). The accuracy in Y p needed to reach these goals ext ends to the 
third decimal place (e.g. IPeimbert et al.ll2003l ; iLuridianal 120031 ; ISteigmanl l2006al ) . Y p deter- 
minations are affecte d by at least thirteen sources of error (see Section [3] and the review by 
Peimbert et al.ll2003l ). In Paper II we obtained a determination of Y p in which most of these 
sources of error were taken into account. In this paper we improve our previous Y p deter- 
mination, considering new atomic data and improved stellar population synthesis models, 
which are likely to further reduce the error affecting our computation. Specifically, four of 
the main sources of e rror are reduced by mean s of the use of: a) the recent He I recombi- 
na tion coefficients by iPorter et al.l (120051 . 120071 ). b) the H I collisional excitation coefficients 
by lAnderson et al.l (120001 . |2002| ) , and c) the correction for underl ying absorption for lines 



redward of 5000 A, based on t he population synth e sis m odels by iGonzalez Delgado et al. 
(120051 ) and the observations by iLeone fc Lanzafamd (119981 ) . 



Y p is determined by means of an extrapolation to Z = of the Y values of a sample 
of objects. Here, the usual normalization X + Y + Z = lis used, where X, Y, and Z are 
the abundances per unit mass of hydrogen, helium, and the rest of the elements respectively. 
he extrapolation is traditionally done in the Y, Z space by assuming a AY/ AZ slope 



(IPeimbert fc Torres-Peimbertlll974l ). More recently it has become common to use AY/ AO, 



where O is the oxygen abundance per unit mass, since the O value is easier to determine. 

In Section [2] we discuss the collisional excitation of the Balmer lines and apply the 
corresponding correction to our H II regions sample. In Section [3] we correct the observed 
line intensities by extinction and underlying absorption, also in this section we determine 
the Y values for each of the objects of the sample, and combining them with a AY/ AO 
relationship we derive the new Y p value. In this section we present the error budget of 
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our determination and a discussion is made where we distinguish between systematic and 
statistical effects. In Section [4] we compare our determination with those by other authors. 
In Section we discuss the cosmological implications of our new Y p value, and in Sections M 
and [7] we present the discussion and our conclusions. A preliminary account of some of the 
results included in this paper was presented elsewhere ( jPeimbert et al.l 120071 ). 



2. Collisional enhancement of hydrogen Balmer lines 

In high-temperature H II regions, the hydrogen Balmer lines can be excited by collisions 
of neutral atoms with free electrons. This effect is generally small, usually contributing less 
than a few percent of the intensity of Ha and H/3 or less, which for most applications can be 
neglected. However, collisional enhancement must be taken into account in the determination 
of Y p for at least two reasons. First, for this task the helium abundance by unit mass, Y, 
in individual H II regions must be known with the highest attainable precision. Second, the 
objects where collisions are more important are those with the highest T e and, hence, the 
lowest metallicities; consequently, they bear a major weight in the extrapolation to Z = 
of the relation between Y and Z. 

In Paper II we used photoionization models calculated with the photoionization code 
Cloudy to estimate the collisional contribution to the Balmer lines in five low-metallicity 
objects. For three objects (Haro 29, 1 Zw 18, and SBS 0335-052 ) we co mputed tailored models 



using version 94.00 of Cloudy, last described by iFerland et al.l (119981 ); for the remaining two 



objectsfNGC 2363 and NGC 346) the mod els were retrieved from previously published works 



(ILuridiana et al.lll999l ; iRelano et al.ll2002l . respectively), and the collisional rate in them was 
estimated from the difference between the total and the case B H/3 emission predicted by 
the models. In the present work, we present improved estimates, which differ from the 
previous ones in several aspects: a) The models for NGC 2363 and NGC 346 have been 
recalculated based on the input parameters described in the original papers, so that the 
collisional contribution could be explicitly computed rather than simply estimated; b) the 
collisional rates have been updated; c) the radiative cascade following collisional excitation is 
now computed self-consistent ly, and d) the model fitting follows a different philosophy with 
respect to Paper II. Points b), c) and d) will be explained in detail in the following sections. A 
further difference with respect to Paper II is that, while the collisional luminosities of Paper 
II models had been computed simultaneously with the models themselves by modifying a 
Cloudy's routine, those of the present paper have been computed by an external program, to 
which Cloudy's output (i.e., the ionization and temperature structures) is fed. The difference 
between the two procedures is, of course, irrelevant from the point of view of results, except 
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for the change in the collisional atomic data used in either case, which will be discussed in 
Section 12.11 



2.1. Updated collisional rates 



I n previous works, it has been estimated that I(H.(3 ) m i/ I(H.a) rn i is about 1/10 ([Davidson fc Kinman 



19851 ; IStasihska &: Izotovl l200ll ; iLuridiana et al.l 120031 ) , and it has been predicted that this 
differential enhancement would increase the observed Balmer decrement, mimicking a larger 
interstellar reddening. The amount of extra reddening due to this mechanism is expressed 
by the collisional reddening coefficient C(Hf3) coL . 



C{H0) 



col 



Log(/(H«) tot //(H/5) tot ) - Log(/(Ha) re 7/(H/?) rcc ) 

-/(Ha) = 
Log((l - J(H/?) col //(H/?) tot )/(l - /(Ha) col //(Ha) tot )) 



-/(Ha) 



Log((l - x (,)/(! - x a )) 
-/(Ha) 



(1) 



where I(X) rec / 'I(\)tot an d x\ = ^(A) co ///(A) to t are the relative contributions of recombi- 
nations and collisions, respectively, to the total emitted intensity in A, and /(Ha) is the 
reddening correction at Ha. Collisional and interstellar reddening add together to yield the 
observed reddening: 

C(R(3) ohs = C(R(3) dust + C(R(3) co \ (2) 

(see Paper II), so that a good estimate of C(Hj3) co1 is needed to properly derive C(Hj3) dust 
from the observed reddening. 



In the present work we have used updated collisional coefficients ( 1 Anderson et al.l 12000 



20021 ) to revise our previous estimates of x a and xp and the derived value of C(H/5) co1 . Two 
major fe atures of the new data deserve to be noted here in comparison to those used in 
Paper II JCarlawavlll99l IVriens & SmeetsHl98oh : 



1. All the collisional coefficients £1(1, ri), to which collisional rates are directly propor- 
tional , are larger than the corresponding values by lCallawayl (119941 ) and lVriens fc Smeets 
( 1l980l ) in the temperature range of interest. 



°The reference to Anderson et al. (2000) and Anderson et al. (2002) contained in Paper II 
as the models had been computed with an earlier version of Cloudy. 



was erroneous, 
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2. The increase of 0(1 , 4) is larger than that of 0(1, 3) by a factor of 2.5 approximately. 



Point 1 above implies larger collisional contributions to all Balmer lines than previously 
estimated. Point 2 implies that J(H/3) co ///(Ha) co / also increases by (roughly) the same factor 
of 2.5 as 0(1, 4)/0(l, 3). This is because, although collisional excitations to any level with 
n > 3 might eventually produce an Ha photon, excitations to n = 3 largely dominate and 
are responsible for more than 80% of the collisional Ha emission in our objects; analogously, 
excitations to n = 4 dominate the collisional H/? emission, yielding more than 90% of it. 
Therefo re, the predicted I(B.P ) m J I (B.a) m i closely follows the 0(1, 4)/0(l, 3) ratio: with the 
data by lAnderson et al.l (120001 . 120021 ) it is now predicted to be larger than 1/4 at T e = 20000 
K and approaching 1/3 at larger temperatures, making it almost indistinguishable from the 
normal Balmer decrement produced by recombinations. As a result, with the updated colli- 
sional coefficients C(Hf3) co1 is predicted to be much smaller than previously estimated, and 
a larger fraction of the observed reddening can be attributed to dust interstellar reddening. 



2.2. Improved radiative cascade 

In addition to using the new collisional coefficients, we have now improved our scheme 
for the calculation of the radiative cascade following a collisional excitation. The collisional 
data used in Paper II were available only as level-averaged coefficients O(n), whereas the 
new coefficients are available as 0(n, /). Therefore, an ad-hoc assumption had to be made 
in Paper II regarding how electrons excited by collisions distribute themselves among the 
sublevels (n, I) of a given level n; this assumption, in turn, determined the radiative spectrum 
following an upward collision. In the present work, in contrast, knowledge of the 0(n, /) 
makes it straightforward to calculate the radiative cascade following upward collisions, and 
no assumptions are required (for low values of n, /-mixing is completely negligible at the 
densities of our objects, N e < 1000 cm -3 ). 



2.3. Model fitting 

Because the collisional enhancement of Balmer lines depends on the interplay between 
the abundance of H° and the temperature structure of the region, good estimates can only be 
given with tailored photoionization models that simultaneously constrain the temperature 
and the ionization structure. In Paper II, we computed several tailored photoionization 
models for each of the three high-T e H II regions SBS 0335-052, I Zw 18 and Haro 29, 
and estimated the collisional enhancement of Ha and H/3 by selecting those models that 



- 6- 



properly fit the observed [O III] temperature, without bothering about the fitting of T [O 
II]. The rationale of that choice was that, because of the Boltzmann factor in the expression 
for the collisional rate, the temperature dependence of collisional rates is very strong: as a 
consequence, we thought that an accurate estimate could only be granted by a good fit of 
the hottest region of the nebula, which in this range of metallicities is the inner one and is 
characterized by the [O III] temperature. 

However, further calculations have shown that, in most cases, the dominant contribu- 
tion to the collisional excitation of hydrogen comes from the outer zones. This happens 
because the outward increase in the fraction of neutral hydrogen outweighs the decrease 
in temperature, exceptions to this behavior will be mentioned in Section [6] and extensively 
discussed in a forthcoming paper (Luridiana, in preparation) . Since the outer zones are best 
characterized by the [O II] temperature, the collisional contribution to Balmer lines is now 
computed based on models that correctly fit the [O II] rather than the [O III] temperature. 
An example of this behavior is given in Fig. [TJ which shows, for one of our model nebulae, 
the relative contribution of each layer to the emission in the [O II] and [O III] lines and 
in H/3coi- The three panels correspond to different integrations, simulating a beam (i.e., a 
point-like slit), a narrow slit, and a slit covering the whole object respectively; depending 
on how the nebula is sampled by the aperture, different layers contribute differently to the 
overall emission. In all cases, we can see that the collisional H/5 emission closely follows the 
emission of the [O II] lines: in other words, H/3 co i preferentially forms in the [O II] rather 
than in the [O III] zone. The effect is more pronounced in the case of a beam, which is a 
good approximation for the data of four of our five objects, but is still seen even in the case 
in which the whole nebula is sampled. 
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Fig. 1. — Normalized contribution per unit radius to the emission in the main [O II] and 
[O III] lines and in H/3 co i in a model for NGC 346, showing that H(3 co i preferentially forms 
in the [O II] zone. The values are normalized to the maximum value in each line. The 
three panels correspond to different integrations, simulating a beam, a narrow slit, and a slit 
covering the whole object respectively. 
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According to this criterion, the fitting philosophy was changed to grant agreement with 
the [O n] rather than the [O ill] temperature. The models of Haro 29 and I Zw 18 used to 
this aim are the corresponding best-fit models of Paper II, w hich already fitte d the o bserved 
T[0 ii] temperatures. For SBS 0335 we used model A of iLuridiana et al.l ( 120031 ). which 
correctly fits the [O n] temperature, instead of the hotter model B, which fits the [O ill] 
temperature and was used in Paper II. For NGC 346 and NGC 2363 (for which in Paper 
II we only gave estimates based on models compiled from the literature) we computed two 
new models with C loudy (version 06.02): the one for NGC 346 is based on model L4 by 



Re 



by 



ano et al.l (120021 ) , and the one for NGC 2363 is based on the best-fit mod el for this region 



Luridiana et al.l (119991 ). In all the models, the updated collisional rates by lAnderson et al. 



( 200ol . 2002 ) were used, and the appropriate corrections for the slit aperture used in the 



observations were made. 



2.4. Resulting estimates of collisional contribution to the Balmer lines 

The models described in Sections 12. 3^ which incorporate the changes described in sec- 
tions [2J] and [221 were used to derive estimates for the collisional contributions. The results 
are listed in Table [lj Based on these values, we recalculated the breakdown of C(B./3) obs in 
terms of collisional and dust reddening for the objects of the sample. The results are listed 
in Table [2J The x a and xp values predicted by our model A for SBS 0335 were not published 
in Paper II. 

The effect of using the new atomic data is an increase in the estimated collisional 
contribution to the Balmer lines. This increase is partially offset by the lower temperature 
of the photoionization models used in this work: some of the models used in Paper II, 
computed to fit the observed T [O ill] values, produced T [O n] values larger by up to 700 
K than the observed values, while the new models correctly fit the T[0 ii] temperature, and 
have correspondingly lower collisional rates. In Table [1] we list the new estimates of x a and 
Xp and the values of T [O ii]. 
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Table 1. Collisional contribution to the total Balmer intensities 21 . 



Paper II This work 



Object 


T e (0 II) 




X/3 




X/3 


NGC 346 


12600 


0.015 


0.004 


0.011 


0.007 


NGC 2363 


13800 


0.030 


0.008 


0.037 


0.027 


Haro 29 


14000 


0.028 


0.007 


0.033 


0.021 


SBS 0335-052 b 


15600 


0.074 


0.021 


0.086 


0.066 


I Zw 18 


15400 


0.060 


0.017 


0.070 


0.053 



a XA = 1 (A) coll 1 (A) text is the ratio between the collisional 
and the total intensit y, calculated wi t h the upda ted colli- 



sional coefficients by lAnderson et a 



photoionization mo dels described in 



2000, 



2002) for the 



Luridiana et al.l (119991 ). 



Relano et all (120021 ) . and Paper II. 



Predictions of model A for the centermost 1.8" x 1" re- 
gion, corresponding to the sum of the three brightest posi- 
tions observed by Izotov et al. (1999). 



Table 2. Observed, collisional, and collision-corrected reddening coefficients^ 







Pap 


er II 


This work 


Object 


C{R(3) obs 




C(K(3) dust 


C(R(3) co1 


C(Rf3) dust 


NGC 346 
NGC 2363 
Haro 29 
SBS 0335-052 
I Zw 18 


0.15 ±0.01 
0.11 ±0.02 
0.00 ±0.08 
0.24 ±0.02 
0.02 ±0.02 


0.02 ±0.02 
0.03 ±0.02 
0.03 ±0.02 
0.07 ±0.05 
0.06 ±0.02 


0.13 ±0.02 
0.08 ±0.02 
-0.03 ±0.08 
0.17 ±0.06 
0.04 ±0.02 


0.00 ±0.00 
0.01 ±0.01 
0.02 ±0.01 
0.03 ±0.01 
0.02 ±0.01 


0.15 ±0.01 
0.10 ±0.02 
-0.01 ±0.08 
0.21 ±0.02 
0.00 ±0.02 



a The uncertainty in C(B./3) co1 was calculated as in Paper II. 
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3. Abundance determinations 

3.1. Adopted line intensities 

In Tables E] and H] we present the observed line intensity ratios of the He I lines relative 
to H(/3), F(A)/F(H/?), for the 5 objects without correction for underlying absorption. 



Table 3. Adopted He I line intensities relative to H/5 a 



NGC 346 b NGC 2363 c Haro 29 c 



1~Tp T linp 

1 _L V J. 1J.AJ.V_j 


F 






F 






J 




F 


I 


3820 









011±0 


001 





014±0.001 





007±0.001 


0.010±0.001 


3889 


0.1820±0.0017 d 


0.0988±0.0019 





174±0 


001 d 





095±0.003 





186±0.001 d 


0.097±0.001 


4026 


0.0171±0.0006 


0.0203±0.0007 





015±0 


001 





020±0.001 





016±0.001 


0.019±0.001 


4387 


0.0045±0.0002 


0.0054±0.0002 





004±0 


001 





005±0.001 





004±0.001 


0.005±0.001 


4471 


0.0370±0.0005 


0.0396±0.0005 





038±0 


001 





042±0.001 





037±0.001 


0.040±0.001 


4922 


0.0100±0.0002 


0.0109±0.0002 





012±0 


001 





013±0.001 





009±0.001 


0.010±0.001 


5876 


0.1143±0.0013 


0.1075±0.0012 





112±0 


001 





110±0.002 





103±0.001 


0.105±0.001 


6678 


0.0336±0.0002 


0.0305±0.0002 





031±0 


001 





030±0.001 





029±0.001 


0.030±0.001 


7065 


0.0243±0.0002 


0.0218±0.0002 





032±0 


001 





030±0.001 





025±0.001 


0.026±0.001 


7281 


0.0073±0.0004 


0.0066±0.0003 





006±0 


001 





006±0.001 





005±0.001 


0.005±0.001 



a The F values and their errors are those presented by the observers. The I values denote the intrinsic 
fluxes after correcting for underlying absorption, collisional contribution to the Balmer lines, and interstellar 
reddening; the errors attached to the I values are only those due to the flux errors, see text. 



3 F values for region A from iPeimbert et al.l (120001 ) . 



C F values from llzotov et al.l ( 19971 ). 
including the contribution due to H8. 



Table 4. Adopted He I line intensities relative to H/5 a 



He I line 


SBS 0335-52 b 


I Zw 18 c 




F 


/ 


F 


/ 


3889 


0.1606±0.0018 d 


0.0948±0.0026 


0.1570±0. 0043 d 


0.0844±0.0052 


4026 


0.0122±0 0005 


0.0179±0.0007 


0.0151±0.0036 


0.0222±0.0053 


4471 


0.0339±0.0007 


0.0401±0.0008 


0.0352±0.0025 


0.0412±0.0029 


4922 


0.0077±0.0004 


0.0094±0.0005 






5876 


0.1168±0.0014 


0.1132±0.0014 


0.0968±0.0028 


0.1016±0.0029 


6678 


0.0322±0.0005 


0.0297±0.0005 


0.0273±0.0019 


0.0290±0.0020 


7065 


0.0453±0.0007 


0.0406±0.0006 


0.0249±0.0016 


0.0263±0.0017 



a The F values and their errors are those presented by the observers. The 
/ values denote the intrinsic fluxes after correcting for underlying absorption, 
collisional contribution to the Balmer lines, and interstellar reddening; the errors 
attached to the I values are only those due to the observed flux errors, see text. 



b F values for the three brightest positions by llzotov et al.l (119991 ): center, 
0".0SW., and 0".6NE. 



C F values for the southeast region by llzotov et al.l (1l999l ). 
including the contribution due to H8. 
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To correct for the stellar underlying absorption in the He I and H I lines we made use of 
synthetic models and observations. The intensity of B./3 was corrected taking into account 
the observed H/3 equivalent widths i n emission, EW Pm .CH.0), and the H/j e quivalent widths, 
EW ab {R(3), based on the models by IGonzalez Delgado et all Jl999l . l2005h (Table EJ). The 
E W n h for the remaining Ba l mer lines an d the He I lines were obtained also from the models 
by IGonzalez Delgado et al.l (119991 . 120051 ) and Cervino (private communication), consistently 
with the adop ted EW^W) values. T h e He I AA 7065 and 7281 A lines are not included in 
the models by IGonzalez Delgado et al. ( 2005 ). therefore we made use of the observations by 
Leone fc Lanzafamel ( 19981 ) to correct for the stellar underlying absorption of these lines. 



After correcting for the underlying absorption and the collisional contribution to the 
Balmer lines presented in Table [TJ we determined the He I intensities (due to recombinations 
and collisions, and affected by optical depth effects) relative to the pure recombination 
H/? line intensities uncorrected for reddening, G(He i)/G(H/3). Once these two corrections 
were made, we combined the C(H/3) dust value with the reddening law by Seaton (1979, 
hereafter S79) and derived the reddening corrected /(He i)/I(H/3) line intensity ratios; note 
that /(H/3) represents the recombination contribution only, while the /(He i) values include 
recombination and collisional contributions as well as optical depth effects. 

In Tables [3] and H] we present the /(He l)//(H/3) line intensity ratios and the errors 
associated with the observational measurements only; specifically, errors associated with 
underlying absorption, collisional contribution to the Balmer lines, and interstellar reddening 
are not included; they will be included in section 3.3 and discussed in section 3.4. We decided 
to follow this procedure because the maximum likelihood method (MLM) used in section 
3.2 to determine the helium physical conditions, requires the errors presented for each line 
to be independent from each other; while the effects of the uncertainties, from any one of 
these three sources (underlying absorption, collisional contribution to the Balmer lines, and 
interstellar reddening), affect the /(He l)//(H/3) ratios in a correlated way. 



For SBS 0335-052 we used the ob s erved F values by llzotov et al.l (119991 ) . We did not 
use the data presented by llzotov et al.l (120061 ) for the brightest 1".56 x 1".04 region for this 
object because, even if the signal to noise is very high, the observations were obtained in seven 
different spectral ranges which did not correspond exactly to the same region of the sky, this 



fact renders these observations useless to determine a very accurate Y value (llzotov et al. 



20061 ). We also did not use the data for the whole object by llzotov et al.l (120061 ) because the 
quality of the line intensity deter minations is considerably lower than the quality of the data 
presented by llzotov et al.l (Il999l ) . 



In Paper I the NGC 346 observations were not corrected for underlying absorption 
in the He I and H I lines because in region A the bright O stars that ionize the H II 
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region were avoided. We consider that this decision was incorrect and that most of the 
continuous emission in region A is due to dust-s cattered light showing the He I and H I lines 



i n abs orption. The reasons are the following: a) lO'Dell &: Hubbard! (119651 ). and lO 'Dell et al. 



(119661 ) showed that most of the non-stellar continuum observed in the visual range of the 
spectrum of H II regions, when the brightest stars are not included in the observing slit, 
is due to dust-scattered light, b) comparing the equivalent width of region A, EW(H/3) = 
250 A, to the e xpected EWCH.0) for a low density pla sma at 12500 K, which amounts to 
~ 1230 A (e.g. lAllerlll984j ; lOsterbrock fc Ferlandl 120061 . and references therein), it is found 
that about 80% of the continuum light is not produced by the nebular gas, c) The color of 
the underlying spectrum obtained by subtracting the nebular continuum from the observed 
spectrum is bluer than that provided by OB stars and therefore cannot be due to A and 
later type stars (the intensity of the observed continuum at differ ent wavelengths can b e 
obtained from the EW of the emission lines presented in Table 5 of lPeimbert et al.l (120001 ). 
and the absolute line intensities presented in Table 2 of the same paper; the intens ity of 
the ne bular continuum at different wavelengths can be obtained from Table 4-9 by lAller 
( 119841 ) and the intensities at different wavelengths f or stars of different spectral types can be 
obtained from Table 8 of the paper by I Co del (119601 ) ) . Since the OB stars were avoided from 
the observations, it can neither be due to them. We conclude from the previous arguments 
that most of the underlying continuum in region A of NGC 346 is due to dust-scattered light 
provided by the brightest OB stars of the cluster. 



3.2. The Y determinations 



The abundance analysis of these objects is based on the combined use of standard 
empirical rel ations and tailored pho toionization models. For NGC 2363, we us ed the models 



descri bed in iLuridiana et al.l (119991 ). while for NGC 346 we used the model by iRelano et al. 



( 120021 ). As mentioned in Section [21 the models for the remaining three objects were presented 
in Paper II. 

In addition to the collisional contribution to the Balmer lines, to obtain He + /H + values 
we need a set of effective recombination coefficients for the helium and hydrogen lines, an 
estimate of the optical depth effects for the He I lines, and the contribution to the He I line 
int ensities due to coll i sional excitation. We used the hydrogen reco mbination coe f ficient s 
by IStorey fc Hummer! (119951 ). the helium reco mbination co e fficien ts by iPorter et al.l (120051 ). 
with the interpolation formul ae provided by IPorter et al.l (120071). and the collisional con - 
tribution to the He I lines by Sawey fc Berrington ( 19931 ) and Kingdon fc Ferland ( 1995 ). 
The optical depth effects in the triplet lines were estimated from the computations by 
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Benjamin. Skillman. fc Smita (120021 ). 



As in Paper I, we used a maximum likelihood method (MLM) based on the He I line in- 
tensities to derive the iV(He + )/N(H + ) values. We produced two sets of models, one where we 
assumed constant temperature (t 2 = 0.000), given by T(4363/5007), and obtained r 3 8 8 9, n e , 
and N (He + ) / N (H + ) , and the other set of models where T(He + ) was an additional variable 
determined also by the MLM, in all cases T(He + ) resulted smaller than T(4363/5007), to 
reconcile both temperatures it is n ecessary to assum e the presence of temperature variations 
characterized by (t 2 ^ 0.000) (see iPeimbertl (119671 ) and Paper I). The resulting values are 
presented in Table where the errors include only those due to: the temperature struc- 
ture, the density structure, the optical depth of the He I triplet lines, and the adopted line 
intensities. 

The MLM solutions with t 2 ^ 0.000 yield lower values of iV(He + )/iV(H + ) than those 
with t 2 = 0.000 (see Table [5]). This change is not due to the variation of the He I recom- 
bination coefficients with temperature, which in this range of temperature is very small. 
The change is mainly due to the higher densities, derived from the MLM when t 2 ^ 0.000, 
which increase the importance of the collisional excitation of the He I lines. The densities 
for t 2 ^ 0.000 of the whole sample are about 67 cm -3 higher than for t 2 = 0.000. The MLM 
densities for NGC 346 and Haro 29 are 4 ± 15 and —7 ± 50 cm -3 respectively, these values 
are too low and support the idea that t 2 has to be larger than 0.000. In Table [5] we have 
adopted for NGC 346 and Haro 29 a density of 10 cm -3 for t 2 = 0.000. The rms density for 
NGC 346 is 14 ± 3 cm -3 (IPeimbert et al.ll2000l ). a reasonable value for this type of objects. 
The rms density is the minimum density that can be associated to an H II region. For real 
nebulae, which always present large density variations, the local density associated with the 
physical mechani sms that produce the line intensities, is alway s considerably higher than the 
rms density (e.g. lOsterbrock fc Flatherlll959l ; IPeimbertl 1 1 9661 ). 



The total helium to hydrogen abundance ratio was derived using the following equation 



iV(He) / N e N(H.e°)dV + J N e N(Re + )dV + J N e N(Re ++ )dV 



N(E) 



ICF{Re 



f N e N(H°)dV + f N e N(R+)dV 
J N e N(Re + )dV + J N e N(Re ++ )dV 
J N e N(R+)dV ' 



(3) 



where ICF(He) is the helium ionization correction factor. The JCF(He) values were ob- 
tained from the Cloudy models and are presented in TableEJ To obtain the iV(He ++ )/iV(H + ) 
values, we used th e 7(4686) //(H/9) value together with the recombination coefficients by 
Brocklehurstl Jl97lh . 



Table 5. Physical Parameters for the H II Regions 





NGC 346 


NGC 2363 


Haro 29 


SBS 0335-052 a 


I Zw 18 




^n 1 1 n 
ZOU ± 1U 


1 27 i 1 n 
lot ± 1U 


oo/i 1 1 n 
ZZ4 ± 1U 


1 «n 1 1 n 

ioy ± iu 


1 rr i 1 n 

loo ± 1U 


EW abs (}I8) 


2.0 ±0.5 


2.0 ±0.5 


2.0 ±0.5 


2.0 ±0.5 


2.9 ±0.5 


/CF(He) 


1.000 ±0.001 


0.993 ±0.001 


0.9955 ±0.001 


0.991 ±0.001 


1.000 ±0.001 


n e [t = (J.UUUj 


1 n _i_ 1 c 

10 ± 15 


one _t_ a t 

205 ± 67 


1 n _i_ en 

1U ± 50 


oon _i_ a 1 


E T _i_ Tn 

57 ± 70 


/ 4-2 n r\f\r\\ 
^3889 = 0.000) 


n 1 n _i_ n on 

0.10 ± U.2U 


1.12 ± U.4U 


1 a e _l n on 

1.65 ± 0.29 


a e _i_ n e 

2.65 ± 0.35 


n n/^ _i_ n n c 

0.06 ± 0.05 


iv (He J / iv (11^ J (t = U.UUU J 


00UO ± 4o 


00 <U ± 1 10 


qr £1 _i_ 1 £;o 

0O0I ± lOz 


Q £ /I O _1_ 1 A P. 

oo4Z ± 140 


o4ol ± 004 


/\7Yu~\ /ATfTJ\f+2 n nnn\b 
1\ {tie) / 1\ {ri){t = U.UUU) 


k 00 _i_ r^n 
o52o ± 5U 


rr A 1 1 ^7/? 

o5o4 ±1/0 


071 n _i_ 1 
iv 1U ± loo 


O 7£ 7 _1_ 1 /I 7 

/ 5 / ± 14 r 


O C 11 _1_ OKA 

oozz ± oo4 


1\ {U)/ iv (11 J (i = U.UUU J 


11 pr 1 1 

11.0 ± l.O 


0.00 ± U.o2 


f .40 ± U. / 


Z.ZO ± U.20 


1 _i_ n 1 v 
1.00 ± U.l I 


0(t 2 = 0.000) c 


139 ± 22 


107 ± 10 


89 ± 9 


27 ± 3 


20 ± 2 


t 2 


0.019 ±0.008 


0.021 ±0.011 


0.024 ±0.008 


0.040 ±0.010 


0.025 ±0.006 


n e (t 2 ^ 0.000) 


58 ±33 


285 ± 92 


61 ±50 


329 ±61 


90 ±80 


73889 (* 2 + 0.000) 


0.10 ±0.20 


1.04 ±0.40 


1.28 ±0.30 


2.56 ±0.35 


0.06 ±0.05 


iV(He + )/A^(H+)(t 2 ^ 0.000) b 


8372 ± 77 


8425 ± 180 


8453 ± 172 


8273 ± 158 


8297 ± 346 


N(Re)/N(R)(t 2 ^ 0.000) b 


8399 ± 79 


8440 ± 181 


8513 ± 173 


8490 ± 159 


8368 ± 346 


N(0)/N(R)(t 2 ^ 0.000) b 


13.5 ±2.0 


10.7 ±0.8 


8.85 ± 1.00 


3.26 ±0.33 


2.16 ±0.25 


0(t 2 ^ 0.000) c 


163 ± 25 


129 ±10 


106 ± 12 


39 ± 4 


26 ± 3 



a Values for the three brightest positions by llzotov et al.l (119991 ) . 
b In units of 10~ 5 . 



c Oxygen abundance by mass, in units of 10 



I 

i— 1 

I 
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From the normalization by unit mass given by X + Y + Z — 1, the iV(He)/iV(H) values, 
the N(0)/N(K) values, and the assumption that the oxygen by mass, O, amounts to 55% 
±10% of the Z value, it is possible to derive Y and O. The N(0)/N(H) and O values are 
presented in Table [SJ while the Y values are presented in Table [6J 

The O/Z = 0.55 ±0.10 value was derived by extrapolating to O = the O/Z yalues for 



the O rion nebula, 30 Doradus in the LMC, and NGC 346 in the SMC derived by iPeimbert 



(12003 ). which amount to (43 ± 5)%, (46 ± 7)%, and (53 ± 8)%, respectively. Note that 
NGC 346 is the O richest object in our sample and that the increase in the O fraction with 
decreasing Z is mainly due to the decrease of the C/O and N/O ratios with decreasing O 
abundance. The error in the O/Z ratio translates into an error slightly smaller than 0.0001 
in the Y p determination. 

The O values in Table are slightly different to those presented in Paper I due to four 
causes: a) the iV(0)/iV(H) values are higher because we took into account the collisional 
contribution to the Balmer line intensities, b) the N(He)/N(H) ratios are higher, c) the t 2 
values are slightly different, because the adopted line intensities are slightly different, and d) 
we assumed that 10% of the O atoms are trapped in dust grains in all objects. 

In Table |6] we explicitly present the AY increase due to the collisional excitation of the 
Balmer lines. Also in this table we present the Y values for t 2 = 0.000 and those derived for 
t 2 7^ 0.000. The Y and Y p values present first the statistical and then the systematic errors 
from all the sources presented in Table [7] (note that the iV(He)/iV(H) ratios presented in 
Table [5] only include a subset of the statistical errors, those due to: temperature structure, 
density structure, optical depth effects and line intensities). For each object we determine a 
Y p value and at the end of the table we present the Y p value for the whole sample, ^(sample). 
The error budget for the Y p (sample) is discussed in the error budget subsection. 



3.3. The Y p determination 

To determine the Y p value from all the objects it is necessary to estimate the fraction of 
helium present in the interstellar medium produced by galactic chemical evolution. We will 
assume that 

AY 

where O is the oxygen abundance by mass. From chemical evolution models of different 
galaxies it is found that AY/ AO depends on the initial mass function, the star formation 
rate, the age, and the O value of the galaxy in question. But AY/ AO is well fitted by a 
constant value for objects with the same IMF, the same age, and an O abundance smaller 
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than~4xl0 3 ( Peimbert et al, l l2007h . Consequently in what follows we will adopt a constant 
value for AY /AO . 



The AY/ AO value derived by iPeimbert et all (120001 ) from observational results and 
models o f chemica l evolu tion of galaxies amounts to 3.5 ± 0.9. More recent results are 



those by iPeimbertl (12003 ), who finds 2.93 ± 0.85 from observations of 30 Dor and NGC 



346, and by llzotov et all (120061 ) who find AY/ AO = 4.3 ± 0.7 from the observations of 
82 H II regions. We have recomputed the Izotov et al. value considering two systematic 
effects not considered by them: the fraction of oxygen trapped in dust grains, which we 
estimate to be 10%, and the increase in the inferred O abundances due to the presence of 
tem perature fluctuation s, which for this type of H II regions we estimate to be about 0.08 
dex (IRelano et al.ll2002l). F r om th ese considerations we obtain AY / AO = 3.2 ±0.7. On the 



other hand IPeimbert et all (120071 ) from chemical evolution models with different histories of 
galactic inflows and outflows for objects with O < 4 x 10~ 3 find that 2.4 < AY/ AO < 4.0. 
From the theoretical and observational results we have adopted a value of AY/ AO = 3.3 ± 
0.7, which we have used with the Y and O determinations from each object to obtain the 
set of Y p determinations presented in Table El 

To determine the Y p average from the whole sample we first need to find the weight that 
should be assigned to each object by considering the confidence we have in each one of the 
determinations. For this we added in quadrature the errors provided by the MLM, presented 
in Table EJ plus all the additional sources of error presented in Table [TJ with the exception 
of the errors associated with the recombination coefficients of both H and He (which will 
affect the sample as a whole and thus will not affect the relative confidence we have in the 
determination from each object). The quadratic addition of these 11 sources of error re- 
sult in: err sll (NGC 346)=0.00311, err sll (NGC 2363)=0.00507, err sll (Haro 29)=0.00475, 
err sll (SBS 0335-052)=0.00583, and err E11 (I Zw 18)=0.00867. And the weights obtained in 
this way, for each object in the sample, amount to: u>(NGC 346)=0.4515, w(NGC 2363)=0.1697, 
w(Haro 29)=0.1927, w(SBS 0335-052)=0.1281, and w(I Zw 18)=0.0579. 

We use these weights to determine the different helium averages (Y) = £\ Y(i)w(i) 
along with the corresponding statistical errors err sta = (Y2i[ err sta('i)u j (i)] 2 ) 1 ^ 2 and system- 
atic errors err sys = Yli err sys(i)w(i)] finally we add the statistical and systematic errors in 
quadrature. 

To compare our results with those of other authors that assume t 2 = 0.000 we have 
computed Y p (t 2 = 0.000) = 0.2523 ± 0.0027, using the five Y and O values for t 2 = 0.000 
presented in Tables and 0, note that this value is not presented in Table Also, from the 
five Y p values for t 2 ^ 0.000 presented in Tabled we derive Y p = 0.2477 ± 0.0029. 



Table 6. Y and Y p values 







\r 
I 


Y 


Y P 




AF(Hc) b 


t 2 = 0.000 


t 2 ^ 0.000 


t 2 ^ 0.000 


NGC 346 


0.0015 ±0.0005 


0.2537 


0.2507 ±0.0027 ±0.0015 


0.2453 ±0.0027 ±0.0019 


NGC 2363 


0.0057 ±0.0016 


0.2551 


0.2518 ±0.0047 ±0.0020 


0.2476 ±0.0047 ±0.0022 


Haro 29 


0.0047 ±0.0013 


0.2577 


0.2535 ±0.0045 ±0.0017 


0.2500 ±0.0045 ±0.0019 


SBS 0335-052 


0.0144 ±0.0038 


0.2594 


0.2533 ±0.0042 ±0.0042 


0.2520 ±0.0042 ±0.0042 


I Zw 18 


0.0114 ±0.0031 


0.2529 


0.2505 ±0.0081 ±0.0033 


0.2498 ±0.0081 ±0.0033 


Sample 


0.0056 ±0.0015 


0.2554 


0.2517 ±0.0018 ±0.0021 


0.2477 ± 0.0018 ±0.0023 c 



a Derived from each object under the assumption that AY/ AO = 3.3 ± 0.7 see text. 

b Increase in the Y abundance due to the collisional contribution to the Balmer line intensities. 

c Equal to 0.2477 ± 0.0029. 
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3.4. The error budget 

Based on the errors presented in Table [5] and Table EJ and the discussion in this subsec- 
tion, we have elaborated the error budget for the whole sample of our Y p determination and 
it is presented in Table [7J The errors in the table are grouped in discrete bins because they 
represent broad estimates. In this table the sources of error are listed in order of importance, 
we will say a few words for some of them. The error budgets of other Y p determinations are 
different to ours for many reasons, they depend on the sample of H II regions and on the 
treatment given by the different groups to each source of error. 

It is difficult to estimate the magnitude of the Y p error caused by each of the thirteen 
sources of error, since there is a significant interdependence between some of them; for 
instance, the uncertainty in the hydrogen collisions will modify the reddening determination 
for each object. To quantify the total effect of each source of error in the Y p determination we 
used the following approach: a) we started with the four errors derived from the MLM (He I 
and H I line intensities, optical depth of the He I triplet lines, collisional excitation of the He I 
lines due to the average density, and temperature structure); b) we ordered the other sources 
of error to be considered one at a time (in order: density structure, helium ICF, underlying 
absorption of the H I and He I lines, reddening correction, collisional excitation of the H I 
lines, O (AY /AO) correction, and recombination coefficients of the H I and He I lines); 
c) we measured how a change in each source affects the determination of all the previous 
quantities; d) we presented the total amount as the estimated contribution to the Y p error 
for this latest source, thus all the cross correlations with the previous sources of error are 
included in the later source. For example, the Y p error produced by the modification of the 
reddening correction due to the hydrogen collisions is presented as part of the error due to the 
hydrogen collisions, and is not presented as part of the error due to the reddening correction; 
also not presented, is the increase in the C(Hj3) formal error due to the uncertainty of the 
hydrogen collisions determinations. 

Estimated in this manner each source of error is independent from the others, even 
the ones labeled as systematic. Consequently the standard deviation of the total error is 
computed by adding in quadrature the standard deviation of all the sources of error. Note 
that we are not assuming that any of the thirteen errors has a normal distribution (even if 
several of the error distributions probably are normal), we are reporting what we estimate to 
be the standard deviation of the Y p error due to each source; also, with thirteen independent 
sources of error included, we expect the final error to be close to a normal distribution. 

In Table [7J we have divided the sources of error in two groups, statistical and systematic. 
The errors labeled as statistical will affect differently the Y p determination for each object; 
therefore, increasing the number of objects in the sample will reduce their final magnitude. 
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On the other hand the errors labeled as systematic, are so in the sense that for any one 
of them our lack of understanding shifts the Y p determination of each object in the same 
direction, for some of them by different amounts, producing an error in the determination 
that cannot be reduced by simply increasing the number of objects; for some of them the 
error can be diminished by selecting a small group of objects where we expect this effect 
to be minimum. The specific equations that should be used to determine the errors from a 
sample, which differentiate the treatment given to the systematic and statistical errors, are 
presented in the fourth paragraph of section 3.3 and were used to estimate the final error of 
the Y p determination. 

As expected the total error obtained from Table [7| is in agreement with the total error 
presented in Table M and implies that the estimates of the different sources of errors in both 
tables are equivalent. 

The most important source of error in the determination of Y p is the collisional ex- 
citation of the Balmer lines. The collisional contribution to the Balmer lines produces a 
non-negligible increase in the Balmer line intensities relative to the case B recombination 
and could introduce a bias in the reddening correction deduced from the observed Ha/H/3 
ratio, both effects affecting the Y determination. For our sample the Balmer collisional con- 
tribution on the reddening correction is practically negligible (see Table [2]). On the other 
hand the collisional contribution to the Balmer line intensities is very important in the Y 
determination (see Table ED- We will give an approximate estimate of the error introduced 
by this effect on our Y p determination. We will consider three sources of error: the collisional 
rates, the radiative cascade, and the model fitting (see section [2]). The published collisional 
rates are expected to become increasingly accurate with time, and the current uncertainty on 
them, estimated from the variation of the Qs published in the last decade, is of the order of 
15% for Q(ls, 3) and Q(ls, 4) in the temperature range of interest to us. These values trans- 
late into almost identical values for I(B.a co \) and J(H/3 co i), which in turn would correspond 
to an uncertainty of about 0.0009 in Y p . The errors on the assumed path for the radiative 
cascade are negligible, since the breakdown of a given Q(ls, n) among sub levels is much less 
uncertain than the absolute value of the Q(ls,n) itself. On the other hand, model fitting 
of the observed H II region is the leading source of uncertainty, because the ionization and 
temperature structures are difficult to constrain to the high level of accuracy that would be 
desirable given the strong dependence of the collisional rates on the electron temperature and 
the neutral hydrogen fraction. In Section [2j we have shown that the collisional contribution 
to the Balmer lines mostly comes from the region where the [O II] lines originate, so that it is 
important for a good model to reproduce the observed T[0 II] temperature. Our estimate of 
the uncertainty is based on the range covered by the collisional contribution in models that 
acceptably fit the T[0 II] temperature. By comparing reasonable photoionization models 
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for the same object we consider that the model fitting introduces an error of about 22% on 
the collisional contribution to the Balmer lines. The combination of the three sources of 
error amounts to about 27% of the collisional contribution to the Balmer lines of our sample, 
which translates into an error of about 0.0015 in Y p . 

We consider the error associated with the collisional excitation of the H I lines as sys- 
tematic because this process has not been studied at length by different groups, hence this 
effect could be systematically lower or higher for all the objects. The error might become sta- 
tistical once the problem is studied further, but in any case we expect it to depend strongly 
on the photoionization models needed for its estimation. 

The second most important source of error is the temperature structure. Most Y de- 
terminations are based on T(4363/5007), but other temperature determinations yield lower 
values, and photoionization models do not predict the high T(4363/5007) values observed. 
These results indicate the presence of temperature variations which should be included in 
the Y determination (see Paper I). The best procedure to take into account the temperature 
structure is to determine T(He + ) based on the maximum likelihood method. The Y abun- 
dances derived from T(He + ) are typically lower by about 0.0040 than those derived from 
T(4363/5007). The error quoted in Table [7] is due to the error in the T(He + ) determina- 
tions obtained with the MLM (see section 3.2). The difference between both temperatures 
is not correlated to the metallicity of the H II region, therefore the systematic error intro- 
duced by the use of T(4363/5007) in the Y determination is similar for objects with different 
metallicities. 

We have estimated the error in the Y p determination due to the adopted effective re- 
combination coefficie nts based on the confidence in the He I line emissivities presented by 



Bauman et al.l (120051 ). The lines used to determine the helium abundances are AA 3820A, 
3889B, 4026AA, 4387A, 4471A, 4921A, 5876A, 6678A, 7065A, and 7281A, where the letter 
indicates the confidence: AA better than 0.1%, A in the 0.1 to 1% range, and B in the 
1 to 5% range. From these values we estima te a systematic err or due to the computed 



emissivities of about 0.0010 in Y. According to lPorter et al.l (120071 ). the error introduced in 
the emissivities by interpolating the equations provided by them in temperature is smaller 
than 0.03%, which translate s into an error in Y p considerably smaller than 0.0001. In our 



preliminary e stima te of Y p (jPeimbert et al.l 120071 1 we used a different interpo lation to the 



Porter et al.l (120051 ) He I atomic recombination coefficients than that given by iPorter et al. 



(120071 ). as well as a slightly different error budget, and obtained that Y p = 0.2474 ± 0.0028. 



The third most important source of error is the extrapolation to zero heavy elements 
content. Fortunately, based on chemical evolution models of galaxies of different types, it 
is found that AY/ AO is practically constant for objects with O < 4 x 10~ 3 , and in good 
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agreement with the observational determinations (see the previous section). One of the 
sources of error in the observational determination is the fraction of oxygen trapped in dust 
grains which has to be taken into account. 



As was demonstrated by lOlive fc Skillmanl (120041 ). the reddening correction is an impor- 
tant source of error. Therefore to estimate the error in the reddening correction we will make 
comparisons among three classic extinction laws and a recent one, th ese laws ar e labeled: S79, 



W58 (Whitford 1958; as parameteriz ed bv Miller fc Mathews! Il972h . CCM89 flCardelli et al. 



19891 ). and B07 (IBlagrave et al.l 120071 ). Fortunately, for our sample the average C(H/3) value 
amounts to 0.09 only, see Tabled For the 4500 to 7300 A region the differences among the 
S79, W58 and CCM89 extinction laws produce differences in the 7(HeA)//(H/3) values of 
our sample of about 0.15%, which correspond to an average difference smaller than 0.0002 
in the final Y p determination. For the 3800 to 4400 A region the differences between S79 and 
W58 produce differences of about 0.4% in the line intensity ratios, which, when combined 
with all the He I lines available, correspond to differences smaller than 0.0004 on Y p . For the 
3800 to 4400 A region the use of the CCM89 law for R v = 3.16, instead of the S79 law used 
by us, gives systematically higher 7(HeA)/I(H/3) values by about 1%, increasing Y p by about 
0.0008 relative to ou r result when all the He I lines presented in Table [3] and Table H] are used. 
Cardelli et al.l (119891 ) made a strong effort to find a simple analytical law to be useful for all 
Rv values, we consider that in the 3800 to 4400 A range the CCM89 law might overestimate 
the extinction because it might n ot be perfectly represe nted by a seventh order polynomial. 
Support for this idea comes from IBlagrave et al.l (120071 ). who find that the CCM89 law for 
Rv = 5.5 overestimates the extinction in the 3030 to 4350A region when compared with the 
B07 law. The B07 extinction law is intended to reproduce the Orion reddening law using 
a modified CCM89 extinction (Rv = 5.5) which only differs from the CCM89 one in the 
3030 to 4350A region. Based on the previous discussion we estimate that the error in the Y p 
determined by us due to the reddening correction amounts to about 0.0007 and we consider 
that most of this error is systematic. Note that for a sample with an average C(H/3) of 0.2 
our estimate of the error in Y p increases to about 0.0016. 

It has often been shown that the correction for underlying absorption in th e H I and He I 



lines has been underestimated (jSkillman et al.l Il998l ; lOlive fc Skillmanl 120041 ; iPorter et al. 



20071 ). We con sider that the best proc e dure t o corr ect for underlying absorption is to use 
the models by iGonzalez Delgado et al.l (119991 . 120051 ). According to these models for those 
objects with EW em (Kf3) > 150A the expected EWab(H[3) amounts to ~ 2A ; for objects with 
EW em (Rf3) < 150A (older objects), the expected EW a b(¥L(3) becomes larger. Therefore 
the correction for underlying absorption for objects with EW em (Hf3) > 150A is inversely 
proportional to EW em (H.j3), while for objects with EW em (H.j3) < 150A the correction, and 
consequently the associated error, increases even faster due to the larger EW a b(H.P) predicted 
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by the models. For samples including a large fraction of objects with EW em (li(3) < 150A we 
ex pect the error for this co n cept t o be larger than ours, moreover to agree with the models 
by I Gonzalez Delgado et al.l (119991 . 120051 ) the average EW a b(H.j3) for any given sample has to 
be equal or larger than 2A . 

In principle, the larger the number of He I lines used to determine a given Y value the 
better. But there are two additional issues that have to be considered for some of the He I 
lines that are not included in our determination: a) the accuracy o f the atomic data for some 
of the lines is lower than that for those lines which we have used ( iBauman et al.ll2005l ). and 
b) the radiative trans fer of the He I singlet lines has to be taken into account for the p — s 
and s — p transitions (IRobbins &: Bernatlll973l . Il974 ). 



It is possible that further work might uncover additional sources of systematical errors 
that would increase the final error. Increasing the sample with objects as well observed as 
those in our sample and with tailor made photoionization models of similar quality as those 
used by us will reduce the statistical errors but not the systematic ones. It is also possible 
that some of the errors that we consider to be dominated by systematics might be reduced by 
further work, for example by choosing a different set of objects where a particular systematic 
effect is expected to be lower, or by providing more detailed observations, or by increasing 
the accuracy of the atomic data determinations, or by providing more realistic models for 
the observed H II regions. 
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Table 7. Error budget in the l^(sample) determination 



Problem 


Estimated error 


Collisional Excitation of the H I Lines 


±0. 


.0015 b 


Temperature Structure 


±0. 


,0010 a 


(AY /AO) Correction 


±0. 


,0010 b 


Recombination Coefficients of the He I Lines 


±0 


,0010 b 


Collisional Excitation of the He I Lines 


±0. 


.0007 a 


Underlying Absorption in the He I Lines 


±0. 


.0007 a 


Reddening correction 


±0. 


,0007 b 


Recombination Coefficients of the H I Lines 


±0. 


.0005 b 


Underlying Absorption in the H I Lines 


±0. 


,0005 a 


Helium Ionization Correction Factor 


±0. 


.0005 a 


Density Structure 


±0. 


.0005 a 


Optical Depth of the He I Triplet Lines 


±0. 


.0005 a 


He I and H I Line Intensities 


±0. 


.0005 a 



Statistical error. 
Systematic error. 
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4. Comparison with other Y p determinations 



The difference between our Y p value of 0.2477 and the 0.2391 value presented in Paper 
II amounts to 0.0086, and is mainly due to: the change in the He I recombination coefficients 
(which produced an increase in Y p of about 0.0040), the change in the H I collisional excitation 
coefficients (which produced an increase in Y p of about 0.0025), the correction for underlying 
absorption in the red He I lines, and the correction of NGC 346 for unde rlying absorp tion. 
In Paper II we used different He I recombination coefficients, those by ISmitsI (119961 ) and 
Benjamin. Skillman. fc Smitd (119991 ). 



It is beyond the scope of this paper to present an error budget for the Y p determinations 
of other authors, but we will discuss some of the reasons responsible for the different Y p 
values derived by other authors. 



There are many differences with respect to the procedure followed by llzotov &: Thuan 

(120041 ). who derived Y p = 0.2421 ± 0.0021. At least two of them are systematic: our use of 



the recent He I recombination coefficients by iPorter et al.l (120051 . 120071 ) , which yield Y values 
about 0.0040 higher than the previous ones, and our use of the recent H I collisional data 
(see Section [2]), which further increase the Y values over the older H I collisional corrections 
by about 0.0025. 



Olive fe Skillman! (120041 ) find a Y p = 0.249 ± 0.009. Our result is in agreement with 
theirs, but our error is smaller. Again there are the systematic differences due to the He I 
recombination data used by both groups and to the estimation of the collisional contribution 
to the H Balmer lines, these two effects probably would increase their result by about 0.006. 



Fukugita fc Kawasakil (120061 ) based on a reanalysis of the llzotov fc Thuan! (120041 ) sample 



of 33 H II regions determined a value of Y p = 0.250 ± 0.004. In addition to a different 
treatment of the underlying H a nd He I absorption the r e are four systematic effects between 
their de termination and ours. iFukugita fc Kawasakil (120061 ) used the He I recombination 
data by iBenjamin et al.l (119991 ). did not take into account the collisional excitation of the 
Balmer lines, adopted the T(4363/5007) value instead of the temperature provided by the 
He I lines to determine Y, and assumed that ICF(He) = 1.000. Including the first two 
effects would increase their Y p value by about 0.009, while the consideration of the third 
effect decreases their determination by about 0.004, the assumption of JC-F(He) = 1.000 
has to be tested with photoionization models, for our sample three of our models showed 
ICF(He) values smaller than 1.000. By using our sample as representative of their sample 
(which might not be true) their Y p value gets reduced by 0.001. Another problem with their 
determination is that it implies a decrease of the He I underlying absorption with metallicity, 
which is not expected; what is expected instead is a decrease of the underlying absorption 
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with an increase of the equivalent width in emission of H/3. 



5. Cosmological implications 



To compare our Y p value with the primordial deuterium abundance D p (usually ex- 
pressed as 10 5 (D/H) P ) and with the WMAP results, we will use the framework of the stan- 
dard big bang nucleosynth esis. The ratio of baryons to photons multiplied by 10 10 , rjxo, is 



given by (jSteigmanl l2006bl ) : 



r] 10 = (273.9 ± 0.3)tt b h 2 



(5) 



where fib is the baryon closure parameter, and h is the Hubble parameter. In the range 
4 < 77 in < 8 (corresponding to 0.2448 < Y p < 0.2512), Y v is related to r] W by (jsteigman 
2006af ): 

Y p = 0.2384 + 77 10 /625. (6) 



In the same r/io range, the primordial deuterium abundance is given by (jSteigmanl l2006al ) : 



10 5 (D/H) P = D f 



(7) 



From our Y p value, the D p value by lO'Meara et al.l (120061 ) . the Q b h 2 value by ISpergel et al. 



( 120061 ). and the previous equations we have produced Table [8J From this table, it follows 
that within the errors the Y p , D p , and WMAP observations are in very good agreement 
with the predicted SBBN values. 
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Table 8. Cosmological predictions based on SBBN and observations. 



Method Y p D p rj 10 Vl h h 2 

Y p 0.2477 ±0.0029 a 2.78±g;|| b 5.813 ± 1.81 b 0.02122 ± 0.00663 b 

D p 0.2476 ± 0.0006 b 2.82 ± 0.28 a 5.764 ±0.360 b 0.02104 ± 0.00132 b 

WMAP 0.2482 ± 0.0004 b 2.57±0.15 b 6.116 ±0.223 b 0.02233 ± 0.00082 a 



a Observed value. 
b Predicted value. 
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6. Discussion 

The effect of collisions on the H Balmer spectrum emerges from this work as the leading 
source of uncertainty for our sample. This uncertainty has in turn three independent sources: 
the theoretical uncertainty on the collisional Qs, the incompleteness of the collisional Os, and 
the uncertainty on the physical conditions of the gas at which collisions occur. 

The last of these factors is probably the most severe. The collisional excitation of the 
Balmer lines is stronger in the hotter zones of the H II regions (which are predicted to be the 
innermost in this metallicity range) and in the less ionized zones (which are the outermost). 
Since the fraction of neutral hydrogen in a nebula varies over a much wider range than the 
Boltzmann factor, which is the main temperature dependence of collisions, the contribution 
of the outer zones dominates, so that T(0 n) is the most appropriate temperature to char- 
acterize collisions when trying to model them. Additionally, in observations covering a large 
fraction of the nebula the emitting volume of the outer parts outweighs that of the inner 
parts, strengthening this conclusion. However, H II regions that are density bounded or have 
a strongly inhomogeneous temperature structure might escape this rule, particularly when 
observed with small apertures. In these cases the uncertainty will be larger. 

It is clear that only detailed observations and tailored modeling will allow us to reduce 
the uncertainty introduced by the collisional excitation of the Balmer lines, particularly in 
the extreme cases; but a better choice might be to avoid critical H II regions altogether, 
by preferentially selecting those H II regions in which collisions are known in advance to 
play only a minor role, i.e. objects with moderate temperatures. Since the temperature 
in H II regions is mostly determined by metallicity, this amounts to saying that metal- 
poor objects are more adequate than extremely metal-poor objects in Y P determinations. 
This conclusion runs counter the common wisdom that the best candidates for primordial 
helium determinations are extremely metal-poor objects (Z < 0.0005), because they permit 
to minimize the error introduced by the extrapolation of the (Y,0) relation to O/H = 0: 
although this is true, a larger uncertainty on Y p is introduced by collisions than the one 
introduced by most of the other sources, including the slope AY/ AO, so the observational 
efforts should be better directed at metal-poor objects. 

Further disadvantages of extremely metal-poor objects in the quest for primordial he- 
lium are that their number is very small and that their H II regions are relatively faint. 
These disadvantages, together with the uncertainty on collisions discussed above, largely 
outweigh the advantage implied by the smaller error introduced by the extrapolation of the 
(Y, O) relation to zero metallicity. Therefore, we are led to the strong conclusion that not so 
extremely metal poor objects, like those in the 0.0005 < Z < 0.001 range, are more appro- 
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priate for the determination of Y p : in this range of metallicity it is possible to find a large set 
of objects with bright H II regions, which will improve the quality and number of emission 
lines available for the determination of physical conditions; furthermore, the temperatures 
of these objects are smaller than those of more metal-poor objects and consequently the 
correction due to the collisional excitation of the H I lines is also smaller. 

This conclusion leads us to another critical point in the appr oach to Y p . It has bee n 
noted that the uncertainty on Y p is dominated by systematic errors (jOlive fc Skillmanl|2004j ). 



As long as this is the case, it is preferable to analyze a few objects in depth and try to correct 
for the systematics than to perform a more shallow analysis of a larger sample, since this 
last method can reduce the statistical uncertainties but not the systematic errors. Hence we 
strongly support the methodological choice of understanding as well as possible the details 
of the objects in a small sample, before directing our efforts towards extending the sample. 
It is only by means of this approach that systematical errors, such as those arising from 
the temperature structure of H II regions, can be gradually understood, corrected for, and 
eventually transformed into statistical uncertainties. 



7. Conclusions 



The new He I atomic recombination coefficients by iPorter et al.l (120051 . 120071 ) increase 
the Y p determination by about 0.0040 with respect to Paper II. 



The H collisional excitation coefficients by I Anderson et al.l (120001 . 120021 ) increase the Y p 
determination by about 0.0025 with respect to Paper II. 

The adoption of temperature variations, instead of the assumption of constant temper- 
ature given by T(4363/5007), reduces the Y p determinations by about 0.003 — 0.006. In this 
paper and in Paper II we did take into account temperature variations and therefore there 
is no systematic difference due to this effect between our two Y p determinations. On the 
other hand, when comparing our results with those of other authors who assume a constant 
temperature given by T (4363/5007), a systematic difference is present that has to be taken 
into account. 

From the analysis of the data studied in this paper we derive a Y p value of 0.2477±0.0029 
for t 2 ^ 0.000, while under the assumption oft 2 = 0.000, we obtain that Y p = 0.2523±0.0027. 
From the Y p given by the WMAP results combined with the adoption of the SBBN, if our 
error budget is correct, our result implies that t 2 has to be larger than zero at more than 
the 1-cr level. 
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It has been argued that the uncertainty on Y p is dominated by systematic errors ([Olive &: Skillman 



20041 ). If this is the case, it is certainly preferable to analyze a few objects in depth and try 
to correct for the systematical errors, than to rely on a statistical analysis of a large sample, 
since this last method can reduce the statistical errors but not the systematic ones. 

From Table [8] it follows that the Y p , D p , and WMAP observations are in very good 
agreement with the SBBN predicted values. Moreover our Y p value, if our error budget is 
correct, provides stronger const raints to some predictions of non-sta ndard Big Bang cosmolo- 
gies than previous studies (e.g. ICyburt et al.ll2005l ; ICoc et al.ll2006l . and references therein). 
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